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(submitted to PRE) 

Patterns forming spontaneously in extended, three-dimensional, dissipative systems are likely to 
excite several homogeneous soft modes (« hydrodynamic modes) of the underlying physical system, 
much more than quasi one- and two-dimensional patterns are. The reason is the lack of damping 
boundaries. This paper compares two analytic techniques to derive the patten dynamics from 
hydrodynamics, which are usually equivalent but lead to different results when applied to multiple 
homogeneous soft modes. Dielectric electroconvection in nematic liquid crystals is introduced as a 
model for three-dimensional pattern formation. The 3D pattern dynamics including soft modes are 
derived. For slabs of large but finite thickness the description is reduced further to a two-dimensional 
one. It is argued that the range of validity of 2D descriptions is limited to a very small region above 
threshold. The transition from 2D to 3D pattern dynamics is discussed. Experimentally testable 
predictions for the stable range of ideal patterns and the electric Nusselt numbers are made. For 
most results analytic approximations in terms of material parameters are given. 
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I. INTRODUCTION 



Spontaneous formation of spatially periodic structures 
on a homogeneous background is ubiquious in nature, 
fascinating to look at, and often hard to understand in 
detail. The periodic structures are almost never ideal. 
Irregularities may be generated in the transient after the 
pattern formation is initiated and anneal after some time. 
Or, in particular in dissipative systems far from equilib- 
rium, they may be the result of an instability of the regu- 
lar, spatially periodic state itself, then often leading to a 
state which exhibits persistent spatio-temporally chaotic 
dynamics. In some systems this is the case arbitrary 
close to the threshold of pattern formation in control pa- 
rameter space. With the help of reduced descriptions 
like phase-diffusion and Ginzburg-Landau like amplitude 
equations, which are to some extent universal (i.e., inde- 
pendent of physical details), several phenomena associ- 
ated with these deviations from the simple periodic struc- 
ture can be explained 

For mostly practical reasons, experimental and theo- 
retical research on pattern formation and dynamics has 
concentrated on quasi one- or two-dimensional systems, 
but most of the results obtained should have a direct 
correspondence also in genuinely three dimensional pat- 
terns. By a genuinely three dimensional pattern (below 
simply 3D pattern) I do here mean a spatially periodic 
structure for which (i) the spatial period(s) are not de- 
termined by the spatial extension of the sample (referred 
to as the class of patterns formed by "competing interac- 
tion" in Ref. 0|) and (ii) for which the spatial extension 
of the sample is in all directions large compared to the 
period(s) of the pattern. This is a stronger conception of 
a "3D pattern" than the one used in Refs. [0j3], which is 



based only on (ii). 

Additional complications in 3D patterns, as compared 
to ID or 2D patterns, arise from the structure and dy- 
namics of defects (dislocations as well as disclinations) , 
which are point like in 2D but line like in 3D. The im- 
plications for dissipative, nonpotential systems, mainly 
those described by the complex Ginzburg-Landau equa- 
tion, have been addressed by several authors But 
there is another particularity of 3D patterns, which has so 
far found little attention: the massive occurrence of ho- 
mogeneous soft modes, which couple to the pattern and 
can drastically change its dynamics. 

By homogeneous soft modes I mean marginally stable 
or slowly decaying homogeneous or long-wavelength per- 
turbations of the homogeneous basic state from which 
the pattern arises. In the abstract sense of the word, 
they are hydrodynamic modes of the basic state. But 
for the sake of clarity the terms "hydrodynamic mode" 
and "hydrodynamics" shall here be reserved for slowly 
relaxing deviations from the thermodynamic equilibrium 
in an unbounded, homogeneous medium, and their dy- 
namics (e.g., the velocity field in a convective flow is a 
hydrodynamic variable). The pattern- forming basic state 
is itself a non-equilibriums state. Thus, although there is 
some correspondence between homogeneous soft modes 
and hydrodynamic modes (in the narrow sense), the no- 
tions are not identical. 

As it has become clear by the investigation of several 
ID and 2D model systems, homogeneous soft modes are 
the key for understanding many of the phenomena oc- 
curring at, or close to, the onset of pattern formation. 
The most prominent example is the mode associated 
with a homogeneous perturbation of the pressure field in 
Rayleigh-Benard convection which leads to a "singular 
mean flow". It is worth noticing that, since it is usu- 
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ally possible to construct self-consistent amplitude equa- 
tions which do not include the effect of homogeneous soft 
modes, their relevance is easily underestimated in the 
theoretical analysis. 

In ID and 2D pattern-forming systems, most hydrody- 
namic modes are damped by the boundaries enclosing the 
system. For example, momentum and heat can usually 
diffuse freely through the boundaries and are stabilized 
by large external reservoirs. Obviously, this mechanism 
is ineffective in systems which are extended in all three 
spatial dimensions. On the other hand, some coupling 
to a reservoir will also be required in 3D in order to sus- 
tain non-equilibrium pattern formation. This could be 
through electromagnetic fields, some matrix embedding 
the active, pattern-forming medium, or some chemical 
reactant provided in excess. But the couplings to the 
reservoirs are highly specific in these cases and stabilize 
only a few hydrodynamic variables. The remaining fields 
do then lead to homogeneous soft modes. As a result, 
several homogeneous soft modes should be considered as 
the rule in 3D, pattern-forming systems. 

For example, when studying the 3D structures formed 
by chemical waves in the Belousov-Zhabotinsky (BZ) re- 
action, the dynamics of the plain BZ reagent does also 
involve convective fluid motion. Since these hydrody- 
namic modes are usually considered to be a nuisance, 
they are suppressed by embedding the reagent in a gel 
[0 . But another homogeneous soft mode excited by the 
pattern, the temperature field (gradients of which are 
probably driving the convection) remains. Since temper- 
ature gradients have a strong influence on the dynamics 
of the pattern , a complete description of the 3D BZ 
reaction should explicitly involve this mode. 

The work presented here is a case study of 3D pattern 
formation in the dielectric regime of electroconvection 
(EC) in nematic liquid crystals. The system was chosen 
because of its easy experimental accessibility. In particu- 
lar, the electric nature of the instability allows to obtain 
patterns with several hundred periods extension in cells 
of a fingernail's size, evolving on the time scale of sec- 
onds. A closely related variant, the conduction regime 
of EC, which always leads to quasi 2D patterns, is cur- 
rently one of the best understood experimental pattern- 
forming systems, on the phenomenological as well as on 
the quantitative level (see the reviews [^2|,^ ) . These ad- 
vantages compensate the inconvenience of dealing with 
rather complicated (electro-)hydrodynamic equations. 

Rather than trying to understand the complex pattern 
dynamics itself, this paper is mainly devoted to the devel- 
opment of consistent reduced descriptions of the dynam- 
ics. Section H sketches the experimental phenomenon 
and the hydrodynamics of dielectric EC, emphasizing its 
3D nature. Approximations used for an analytic or semi- 
analytic description of dielectric EC are introduced in 
Section III, thereby discussing the linear stability prob- 
lem. In Section |l^ the 3D amplitude formalism for di- 
electric EC is derived. Close to the threshold of EC 
and in a liquid-crystal slab of large but finite thickness. 



the pattern dynamics becomes essentially 2D. The cor- 
responding equations of motion are derived from the 3D 
formalism in Section In Section VI the stabilit y of 
ideal periodic patterns is investigated and in Section VII 
a general scenario for the transition from the onset of 
dielectric EC to fully 3D pattern dynamics with increas- 
ing external stress is developed. Section VIIl discusses 



possible experiments based on electric Nusselt number 
measurements and Section [X summarizes the results. 
Appendix ^ contains some analytic and numerical re- 
sults for coupling coefficients. Appendix ^ compares two 
different methods for integrating multiple, homogeneous 
soft modes into the amplitude formalism in a general 
framework; one method is used in the main text. 



II. SOME PHENOMENOLOGY OF 
ELECTROCONVECTION 

Notice that below some points are oversimplified in 
order to ease intuition. For comprehensive reviews of 
EC see Refs. [p"2|- p{, f or introductions into nemato- 
hydrodynamics Refs. |]l5| , |l6| . 



A. Basic phenomena 

In the typical experiment a nematic liquid crystal with 
negative dielectric anisotropy is sandwiched between a 
pair of transparent, parallel electrodes (separation d ^ 
20 — 50^m, area ~ Icm^). By a special treatment of 
the electrode surfaces, the nematic director n (the lo- 
cally averaged molecular orientation; |n| = 1) is forced 
to align parallel to the electrodes in some preferred direc- 
tion which shall here be identified with the a;-direction {z 
be normal to the electrodes, y normal to x and z). An 
ac voltage E^zd cosojt is applied at the electrodes. In 
the conduction regime at frequencies below the cut-off 
frequency tOc, the first instability to be observed as the 
voltage is increased is towards a pattern of convection 
rolls called Williams domains . 

At higher frequencies a different kind of structure peri- 
odic along X is found. Compared to Williams domains it 
has shorter wavelength and decays faster after switching 
of the voltage (fast turnoff mode). At least two concur- 
ring mechanisms have been proposed for this high fre- 
quency mode: the dielectric EC ||l8|] , which depends es- 
sentially on the anisotropy of the nematic (its threshold 
diverges at the nematic-isotropic phase transition |p^), 
and the isotropic mechanism pO| , pl| where the liquid 
crystal's anisotropy is not essential for the convection 
mechanism itself but only for selecting a preferred modu- 
lation direction. It has a finite threshold at the nematic- 
isotropic phase transition as its characteristic signature 
p2| . The two linear modes have the same symmetry 
and do in principle mix, but generally the corresponding 
thresholds can be assumed to be sufficiently separated to 
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consider the mechanisms isolatedly. The isotropic mode 
is thought to be located mainly near the electrodes, while 
the dielectric mode is maximal at mid plane. Unfortu- 
nately, it is not always clear which mode is actually ob- 
served. At least in some cases the dielectric mode could 
be identified by the good match of the threshold curve 
with theoretical predictions (e.g. [^). The isotropic 
mechanisms will not be considered here. 

For voltages slightly higher than the threshold of di- 
electric EC, the formation of the chevron superstructure 
is observed: defects (dislocations) in the pattern of con- 
vection rolls accumulate along lines oriented in y direc- 
tion, such that the topological charge of the defects alter- 
nates from line to line. Between the lines, the convection 
rolls are rotated and the nematic director is twisted, al- 
ternately clock- and counterclockwise [Q. The observa- 
tion of chevron patterns in the conduction regime of EC 
with homeotropic director alignment p5| , |26t shows that 
this scenario is not restricted to a particular convection 
mechanism. 



B. Hydrodynamic equations and material 
parameters 

EC in both the conduction and the dielectric regime 
result from the interaction of electric field, space charges, 
mass flow, and the nematic director via the Carr-Helfrich 
mechanism: Spatial modulations of the director ori- 
entation are amplified by an inhomogeneous mass flow 
generated by electric volume forces on space charges 
which accumulate due to inhomogeneous electric currents 
in the inhomogeneous director field. Thus Maxwell's 
equations ^ (in the quasi-static approximation ^ curli? = 
curli? = 0) and the balance equations for charge, mo- 
mentum, mass (continuity equation), and the torque act- 
ing on n have to be taken into account. They contain 
several material parameters: the conductivities a\\ ^ 
a± = O{10~^ ■ ■ ■ 10~^ m~^) for electric currents par- 
allel (II) and perpendicular (_L) to n respectively (they 
vary on a large range depending on purity and doping, 
while (7||/crj^ changes only little), the dielectric constants 
^ e||,e_L — 0(eo) (the quantities aa (T[| — (t± = 0{(t±), 
Ea := e|| — e± = 0{e±) measure their anisotropics), the 
fiexoelectric constants 61,63 = O(10~^^ — 10~^^Cm~^) 
(6+ := 61 -|- 62, 6_ := 61 — 62), the diffusion constants for 
(ionic) charge carriers 0(10^^^ m^ s~^) =: Dp, the mass 
density pm = 0(10"^ kgm""^), the five independent vis- 
cosities ai, a<^ = 0(0.1 Nm~^ s) (ag = a2 + 0^3 + as, 
7l = "3 - a2, 72 = "3 + ^2, 2771 = -a2 -I- Q!4 -I- as, 

112 — I2 + Emd the curvature elasticities of the direc- 



tor field ^22 < fell < = 0(10"" N). I also include 
the "dynamic fiexoelectric effect", which was predicted 
p8| , p9[ on the basis of a systematic rederivation of ne- 
matohydrodynamics, but, to the authors knowledge, has 
not been detected, yet. It is characterized by a param- 
eter and leads to additional dissipative contributions 
in the charge, momentum, and torque balance equations. 



C. Dimensional analysis 



With the exception of the charge relaxation time tq :— 
e±/a± = O{10~^ — 10~^s) the nematohydrodynamic 
equations (without external fields), being derived as a 
limit of large time and length scales (though typically 
valid down to molecular scales), do not set any time or 
length scale by themselves. Instead, one finds basically 
three types of diffusivities: for charge (Dp), director ori- 
entation [e.g. -Dd.stat = ^33/71 = 0(10"^" m^s"^) for 
static, Dd^dyn = (fc33?7i)/(7i?yi - ai) = 0(10"'' m^ s^^) 
for dynamic deformations; notice that (static) flexoelec- 
tric effects do not introduce a new diffusive scale since 
6^^3/eo ^ fcss], and momentum [e.g. Dp — (71771 — 
aD/^-jiPm) = O(10~^m^s~^) along n\. The over- 
damped limit p„i 0, Dp — > cx) is generally a good 
approximation. Charge diffusion is not essential for the 
Carr-Helfrich mechanism and is typically screened out. 
Then orientational diffusion sets the only diffusive scale. 

With an externally generated electric ac field 
Eqz cos ojt two additional time scales are introduced: The 
period 27tuj~^ and the "director time" = 7i/(e^i?g) 
[the more intuitive choice Td :— 7i/(eai?o) would suggest 
that ea = is singular for EC, which is, for the convective 
modes themselves, not the case]. 

In the conduction regime the charge densities oscillate 
with the frequency of the applied field, while the director 
orientation is mostly constant. This imposes a condition 

Td > iU-' > To (1) 

on the three time scales. The finite sample thickness d 
determines the wavelength A of the convection pattern 
and leads through a condition d ~ A ~ (fd.stat TdY^^ to 
a voltage threshold = d^E^ ^ fc33/e_L for the onset 
of convection (see Ref. [Q for a good analytic formula). 
A lower limit for the sample thickness is given via rela- 
tion (P by (P > Dd,statTQ- At frequencies higher then the 
cutoff frequency cUc ^ Tq^ the conduction mechanism is 
also disabled or at least supersede by the dielectric mode. 



iln MKSA units. 

^We will only allow for homogeneous magnetic fields. 
^For convenience, the factor eo is absorbed into ey and ex. 
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In the dielectric regime director and fluid flow oscillate 
with to, which leads to a condition 



Td 



< 



(2) 



The charge distribution is at high enough frequencies 
{uj~^ ^ To) mostly constant in time. This is not ac- 
tually necessary for the dielectric mechanism to be ef- 
fective [Q, but typically dielectric EC is supersede by 
the conductive mode at lower uj. The threshold for the 
onset of EC is now given by a condition Td ^ uj^^ (or 
E'^e±/ji ^ uj), i.e., the lowest Eq compatible with re- 
lation (||). The wavelength A of the critical mode can 
under some conditions be ^ d [ pO| , but for typical mate- 
rials used it is A ~ (D^.tiyn/w)^ , at least as long as this 
length is smaller than d and larger than the Debye screen- 
ing length {DpToy/^, i.e., utq < Dd,dyn/Dp = 0(10^), 
where charge diffusion becomes important. 

Thus, the length scales given by the spatial period of 
the pattern A and the sample thickness d are usually 
independent and easily separated (A d) in the dielec- 
tric regime, either by increasing d or by simultaneously 
increasing to and the conductivities, while leaving the 
secondary control parameter ujtq constant. Since strong 
doping may affect the nematic material parameters and 
the nematodynamics at high frequencies is not fully un- 
derstood, the program carried out below is best seen as 
the limit of thick cells. The theory should accurately de- 
scribe typical experiments in cells with d > 10 A. To ob- 
serve fully three dimensional patterns, thicker cells might 
be required (see also Section 



VII) 



the distance from the threshold of pattern formation in 
the control-parameter space of the underlying system. 

It will be shown in Sec. |^ that, as a direct consequence 
of the separation of length scales in the dielectric regime, 
the range of validity of Eq. (||) is highly restricted. Al- 
ready for values of e' of the order A'*/c?'* corrections to 
Eq. (^) must be taken into account. For e' ~ A^/c?^ the 
2D description breaks down completely. But the con- 
vective dynamics can then still be described in terms of 
the 3D modulations of the complex pattern amplitude 
A{x,y, z,t) (defined by Eq. (^ below), which is coupled 
nonlinearly to several homogeneous soft modes. It is 
therefore natural to derive first the 3D amplitude dynam- 
ics which can then be reduced further to a 2D description 
in a subsequent step. 

Julien, Knobloch, and Tobias were the first to im- 
plement the idea of deriving a reduced description for the 
z-dependence of the amplitude of patterns with X/d ^ 1 
as an intermediate step in the theory, and also the first to 
observe that this method significantly eases the restric- 
tion of the control parameter to values close to thresh- 
old. Their calculation do, however, not involve in-plane 
modulations of the pattern and the resulting excitation 
of homogeneous soft modes. Several results concerning 
the 3D description of dielectric EC and its reduction to 
2D are derived in an unpublished work by Lindner [ p3[ , 
which is quoted here whenever necessary. 



B. Linear stability in 3D 



III. APPROXIMATION METHODS AND LINEAR 
THEORY 

A. 2D vs. 3D amplitude formalism 

Below we will develop the amplitude formalism for the 
pattern dynamics in the dielectric regime, i.e., obtain the 
laws of motion of amplitude and phase of the spatial mod- 
ulations as described by the complex pattern amplitudes 
A'{x,y,t) or A{x,y,z,t), respectively. 

The basic state in the experimental cell is anisotropic 
and inversion symmetric and the primary bifurcation is 
supercritical (forward) towards a steady-state pattern 
with a single critical wave vector. Hence, the most el- 
ementary description of the pattern dynamics is give by 
the time-dependent Ginzburg-Landau equation in 2D, 



rd,A'^(e' + edl + eyd'-9'\A'n A' 



(3) 



important physical properties of which are reviewed in 
Refs. ||3l| , ^ . The real, positive coefficients r, ^x, £,y, and 
g' have magnitudes corresponding to natural scales of the 
system (e.g. the pattern wavelength for ^x, Cy) and can 
be calculated from the underlying hydrodynamic equa- 
tions. The small, dimensionless parameter e' measures 



The starting point for setting up the 3D amplitude 
equations is to calculate the linear threshold Ec, critical 
wave number Qc and critical eigenvector (i.e., the I-k juj- 
periodic time dependence of the hydrodynamic fields at 
threshold) ignoring any spatial variations along z. Sev- 
eral linear stability calculations of this type have been 
carried out [|l|,|o|j34|jl|] . 

In experiments, the critical wave vector is always found 
to be parallel to the orientation of the nematic director 
in the basic state (the x direction). The linear problem is 
thus effectively one-dimensional, with trivial, sinusoidal 
variations along the remaining x direction, and is much 
easier to solve than the 2D problem including variations 
and boundary conditions along z. This technical advan- 
tage of the 3D approach, which is of course not restricted 
to EC, remains effective also in the subsequent calcu- 
lations of the coupling coefficients in the 3D amplitude 
equation. 

For analytic as well as numerical calculations it is con- 
venient to use a truncated Fourier expansion of the time 
dependences of the hydrodynamic fields, assuming them 
to be 2^1^} periodic. In the simplest cases truncated 
at lowest orde r (i.e.. including constant and sin/cos 
contributions) |L3| , p3[ or including the sin/cos 2<x't modu- 
lation of the induced electric potential in order to bet- 
ter model the interplay between electric charges and 
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fields. With these truncations, the stability problem for 
sinusoidal excitation can be solved explicitly (see Ap- 
pendix Notice, however, that this "lowest order" 
or, respectively, "second lowest order Fourier approxima- 
tion" involves some arbitrariness in the choice of variables 
and does not correspond to any physical limit. Numerical 
convergence (5% accuracy) requires inclusion of at least 
the third harmonic. The actual time dependence of, e.g., 
the director field, depends on uj and is non-trivial even 
as a; ^ oo Q . 

Some results presented here rely on the first or second 
lowest order Fourier expansion of the induced electric po- 
tential <&, the nematic director expressed by and Lp 
such that n = UzZ + (1 — nlY^'^c, c := (cos 1^9, sin 0), 
the velocity field v, and the pressure P. The other 
hydrodynamic fields are treated implicitly. For the 
representation of the dielectric mode itself, v is ex- 
pressed in the divergence-free form v = {—dy,dx,0)g + 
{dxdzjdydz, —d1 — dy) f and the pressure is ehminated. 

As a natural consequence of Td ^ the relative 

phases of electric, director and velocity fields in the lin- 
ear eigenvector are shifted by angles 0(1). Remarkably, 
the phase shift between the lowest Fourier mode of di- 
rector oscillations and external field is 7r/4, for the first 

f3[ and second-lowest-order Fourier approximation (see 
ppendix ^) exactly and only slightly perturbed (< 1% 
for lutq > 2) when higher Fourier modes are included. 
No simple physical explanation for this result should be 
expected, since it holds only at the critical (most unsta- 
ble) wavenumber. Experimental observations | |35[ | seem 
to agree with a value 7r/4 for the phase shift even better 
than the comparison with Galerkin calculations includ- 
ing z-dependence js^j, which had been carried out as a 
test of the dielectric model of EC. 

We define the pattern amplitude A such that the am- 
plitude riz.c of the cosut Fourier component of the direc- 
tor tilt oscillations (which is in phase with the applied 
voltage) has a spatial dependence 

nz,c = M^i Hi t) ex.p{iqcx) + c.c . (4) 

Assuming as usual e, dx, dy, dz to be small and dis- 
carding contributions beyond the lowest nontrivial order, 
the linear part of the 3D amplitude equation assumes the 
form 

TdtA = (e + ildl -f ildl + ezdl) A. (5) 

As conventional, e := {E^ — E^)/E^. The coherence 
lengths S,y, S,z turn out to be A and the relaxation 
time T is of the order of the charge relaxation time tq js^] 
(s. Appendix ^). 

The horizontal boundary conditions for A are simply 

^ ^ at z = ±d/2. (6) 

Contributions from derivatives of A and nonlinear con- 
tributions to the boundary conditions are of higher order 
and can be discarded. In particular, as is well known, the 



distinction between free and no-slip boundary conditions 
for the velocity field plays no role at this point. With the 
realistic no-slip boundary conditions for the velocities, 
the relative magnitude of the hydrodynamic fields in the 
2D linear eigenvector (including x and z variations) lo- 
cally deviates from the ID eigenvector (only x variation) 
only in a boundary layer of thickness ~ A, an example 
of which is shown in Fig. |l|. This boundary layer might 
provide a problem for numerical approaches directly us- 
ing 2D eigenvectors in thick cells with no-slip boundaries, 
notably when Galerkin approximations are used. 



C. Linear stability in a cell of finite thickness 



Assuming as usual the lateral {x, y) extensions of the 
cell to be large compared to its thickness d, the trivial 
solution A = of Eq. (||) with boundary conditions (||) 
becomes unstable at e = := (tt^z/c^)^ [i-e., e' = e — Cd in 
Eq. (^)] with a critical mode A ~ coa{nz/d). The small 
threshold shift Ec — > {l+ed/2)Ec due to the z-variation is 
rather uninteresting by itself (there is also a shift A'^/d^ 
in qc by a discarded contribution to Eq. (|^) of the form 
idxdzA), but the effect provides, for example, a simple in- 
terpretation of the small gap {AEq — 3edEc/2) between 
the critical mode and the lowest z-antisymmetric mode 
(i.e. A ~ sm{2'Kz/d)) reported in pTf . 



D. Flexoelectric effects 



A short remark about flexoelectric effects, which are 
generally difficult to isolate experimentally, is in place 
at this point. The high symmetry of the linear problem 
in 3D does not allow flexoelectric effects: Ec, qc, t, and 
S^x are independent of the flexoelectric coefficients. Most 
of the remaining linear and nonlinear coefficients con- 
tain flexoelectric contributions. With our choice of vari- 
ables these contributions are, except for some dynamic- 
fiexoelectric terms, all indirect: fiexoelectric effects ex- 
cite additional, "slaved" contributions in the subspace 
orthogonal to the critical eigenvector, which, again by 
fiexoelectric effects, feed back into the dynamics of the 
amplitude of the eigenvector. In contrast, the contri- 
butions not depending on flexoelectric coefficients are, 
except for ^xi all direct: no excitation of slaved degrees 
of freedom is involved. Therefore flexoelectric effects are 
separated in a natural way from the standard dynamics. 
This might provide methods for measuring the flexoelec- 
tric coefficients in a way not sensitive to parasitic bound- 
ary effects. Previous calculations involving flexoelectric 
effects [^,Q|9| where restricted to the conventional, 
"static" flexoelectric contributions and concentrated on 
the determination of critical mode and voltage. 
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IV. THE 3D AMPLITUDE EQUATIONS 
A. Method 

Before discussing the homogeneous soft modes relevant 
for dielectric EC, some comments on methodology are re- 
quired. 

Multiple homogeneous soft modes excited by a pat- 
terning mode have, to my knowledge, first been intro- 
duced by Plant and Pcsch Q. But their description 
requires the soft-mode amplitudes to be constant along 
all but one spatial direction. This limitation seems to be 
partly due to the procedure by which the equations were 
derived. There are two popular philosophies for this pro- 
cedure (see below and Appendix ^ , which shall here be 
labeled as "order parameter" method and "center mani- 
fold" method. The two methods usually (for at most one 
homogeneous soft mode) differ only in the way in which 
the problem is formulated and solved, but lead to the 
same results. The association of existing general prescrip- 
tions for deriving amplitude equations with the former 
(e.g. ||l|-||,||) and the latter philosophy (e.g. ||,|6|, 
the Chapman-Enskog approach [ ^J48| for the derivation 
of hydrodynamics from statistical mechanics is also of 
this type, see RefsJ M 50j ) is therefore not always conclu- 
sive (see also |l|,|l||2|nL When using the "order param- 
eter" method to obtain reduced equations, each spatial 
Fourier mode of the physical state is projected onto the 
(adjoint) slowly decaying linear eigenmodes of the basic 
state with the corresponding wave vector. Using the "cen- 
ter manifold" method, the projection is always onto the 
(adjoint) eigenvectors for (typically neutrally stable) ho- 
mogeneous perturbations. When there are multiple slow 
modes at a single wave vector (usually q = 0), the result- 
ing reduced equations differ, as is shown in Appendix 
Plant and Pesch Q seem to be using the "order param- 
eter" method, which leads to problems in more than one 
spatial dimension. Here, the "center manifold" method is 
used to derive the nonlinear extensions of the amplitude 
equation (^ including homogeneous soft modes. 



B. Derivation of the soft mode equations 

Some particularities of the problem under considera- 
tion have to be taken into account: The quasi-static ap- 
proximation of electrodynamics curl E — — curl B and 
the approximation of an incompressible fluid ViJ — 
both lead to additional homogeneous "soft modes" : Since 
no time derivatives of electric potential (electric field) or 
pressure occur in the basic equations as they are used 
here, all their temporal Fourier modes are in the kernel 
of the linear operator L(0, 0) (s. Appendix^). When the 
viscid limit pm —>■ becomes effective, the same applies 
for the oscillating part of the velocity field. This is the 
case when spatial variations occur on scales smaller than 
(a;L»p)i/2 = {Dp/DaynY^^X = O(105/2)A. Formally we 



shall assume length scales to be larger than this. But, 
since, at least in this case, velocity and pressure oscilla- 
tions do not feed back into the remaining dynamics at 
lowest order in the derivatives, i.e., there are no contri- 
butions 0(/9,~^), the description should be good also on 
smaller scales. 

To simplify the problem further, the equations for the 
soft modes are here calculated only in the lowest order 
Fourier approximation. The "slaved" modes being elimi- 
nated are then the slowly varying average director tilt , 
which is stabilized by the applied electric field through 
the dielectric anisotropy ea (assumed to be negative here- 
after), and oscillations of the director, which are viscously 
damped. Finally, in anticipation of corresponding bound- 
ary conditions, only small deviations from the basic state 
v,ip,nz,^ — shall be considered for now. 

In order to obtain a consistent truncation of the soft- 
mode equations, recall that the only non-diffusive scale 
in the hydrodynamic equations is the charge relaxation 
time tq. Assume lotq to be fixed. This also determines 
Ec ~ uj^^'^ for given tq and we will assume Eq ~ Ec- 
The elimination the fast modes (n^ and oscillations of 
n) becomes more efficient when uj and Eq increase or, re- 
spectively. To decreases. Then, in the limit of small tq (or 
large (t\\ , with fixed cry/cr^), the time scale tq drops out of 
the equations. A purely diffusive scaling for the deriva- 
tives {d1 ^ dy ^ d"^ ^ dt) is retained, without making 
any a priori assumptions about the actual scaling laws of 
typical lengths and times, which may be different. This 
approximation breaks down when length scales become 
shorter than (L'd.statTo)^^^ or 

It turns out that with these approximations the only 
relevant modes in the fast subspace (i.e. R in Ap- 



pendix B 3 ) are Uz and its temporal modulations as given 

by 

nz{x,y,z,t) = 

{k22 - kii) dyUy + asVx 



2dz 



E^ea 

egEg^r - (e+ - 27iC^) g^^p + ol^Vz 
Elta 



(7) 



26 (ca^x^o - e^dyUy) 



s;3f2 



(i^gCa cos(wi) — AjiU) sm{ujt)) . 

The electric potential has been decomposed as $ = 
$0 -I- 2^rCOfiujt — 2^iSmu!t. The parameter b :— 
-^0 ^a/(3 ^-o + 16 7^w^) measures the strength of the 
excitation of the oscillatory part of n^. It is nu- 
merically small [in the standard material MBBA (p- 
methoxybenzilidene-p'-n-butylaniline) 6 ~ 0.01]. The re- 
sulting description for the soft mode dynamics is given 
by the Eq. (|40|) on page |l^ (for the terms containing A, 
see Section |IVF| below) . 
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C. Comments on the soft mode equations 

Most of the terms in Eq. ( ^ ) reproduce linearized ne- 
matohydrodynamics. Equations ( 40a -BOc ) derive from 
the charge balance equation [Eqs.(40t ,40c) have been 
multiplied with Eq], Eq. (40d]from the angular momen- 
tum balance on ft, Eqs . (40c-40g) from the Navier-Stokes 
equation and Eq. ( 40h ) is the unchanged continuity equa- 
tion. In Eq. (40d) a term XaHy(p has been included, 
which describes the action of a magnetic field in y di- 
rection. It will be used in Section The contributions 
resulting from the elimination of fast modes are under- 
lined. 

Remarkably, the equations are mostly independent of 
the strength of the external field (a factor Eg can be ab- 
sorbed into the definitions of $r and <I>i), although it is 
the cause for the excitations of the slaved modes. As 
an example, consider the mechanis m for the reduction of 
viscosity by the term a2d'^Vz in Eq. (40g) (recall a2 < 0): 
Shear forces ~ OL2dxVz excite n^. This leads to polariza- 
tion charges ~ dxEoeanz. The electric field ~ Eq acting 
on these charges generates bulk forces on the fluid. On 
the other hand, the excitation of rt^ is damped by electric 
forces ~ SaEg and the factor Eq cancels out. 

The soft mode equations reflect the non-equilibrium 
character of the basic state. For example, if Onsager 's re - 
lations would hold, th e coe fficients of dxdzVx in Eq. ( |40^ ) 
and of dxdzVz in Eq. ( [40e| ) would be the same. At low lu 
the basic state is even unstable. The mechanism corre- 
sponds to EC in the conduction regime. Ignoring flexo- 
electric effects, the cutoff frequency lUc above which the 
basic state stabilizes is given by 



("2 e± CT|| - a2 e\\ cr± - Ca Vi 



(8) 



which reproduces the result of direct stability calcula- 
tions using (effectively) the same Fourier truncation . 
Only the threshold field for the Williams domains is too 
small to be resolved by Eq. (EO). 



D. Nonlinear extensions 

Of course, constant values can always be added to 
any of the soft modes by a Galilei transformation, a 
rotation, or a gauge transformation. The problem of 
adding nonlinear contributions (e.g. advection terms) to 
Eqs. (Eg) such that they become formally invariant un- 
der these transformation is easily solved. The solution 
is not unique, but it can be seen by inspection that the 
precise form of the nonlinearities does not matter under 
the following conditions: 

1. Dynamics is such that, in fact, diffusive scaling 
holds. In particular this implies that, if L ^ \ 
is the typical length scale (i.e. dx,dy,dz ~ L^^), 
the variations of v, ^r, ^i, f/L, and ^q/L over L 
scale in the same manner as L ^ cx). 



2. Variations in if over L are much smaller than one. 
It is not necessary to specify the scaling relation of 
ip and L. When L is determined through the bulk 
dynamics, ip may actually vary by 0(1) over the 
sample. 

It should be noticed that, although the second condition 
is satisfied for many problems of pattern dynamics, it is 
too strong for disclinations (line defects) in the director 
field at any distance R from the core of the disclination: 
On the typical length scale L — R variations of (p are 
0(1). Then, for example, nontrivial nonlinear contribu- 
tions from the velocity field v , like a term of the 
form {dxVz){dy(p) in Eq. (40g), might have to be included. 

For the part describing the curvature elasticity in 
Eq. ( [40d| ) , the fully nonlinear corrections in ip have been 
calculated by extending the "center manifold" method 
to nonlinear contributions and requiring rotation invari- 
ance. The result has the same form as close to equilib- 



-fidtp = h3c± ■ {dl + dy)c 

+ {k[^ - fc33)(c± • V)(V • C) + k22Cl_ -dlc^... (9) 

[c^ := (— sin cos (p, 0)], however, with :— k\\ -|- 
2be^_/ea. Although an additional term proportional to 
(V • c^)(V • c) would be thinkable, it does not occur in 
the present approximation. For a rotationally invariant 
description of the dynamics of the convection pattern, 
it is useful to go over to a representation in terms of 
A{r) ~ A{r) ex.p{iqcx) as in Ref. [5^. 



E. Stability of the twist mode with respect to x 
modulations 



Another point to notice is that below the threshold of 
dielectric EC (i.e. with ^ = 0) the system can, in our 
approximation, never destabilize in such a way that p is 
excited but dycp = 0, even when allowing for flexoelectric 
effects and arbitrary z dependencies. Without y mod- 
ulations, the in-plane director couples only to Vx, and 
this interaction is not affected by the elimination of the 
fast modes and hence relaxational. This is remarkable 
because such modulation instabilities of p below the EC 
threshold are apparently observed in experiments [|l9|j54[ 
(the "inertial mode" [ p5[ which was proposed as an ex- 
planation has the wrong symmetry). A far-fetched but 
possible explanation would be that the modulations in 
p> are generated through a mechanism which is similar 
to the one that generates the chevron superstructure in 
the dielectric regime, however, invoked by the convection 
rolls of "isotropic" EC (see Section II A). The isotropic 
mechanism is expected to become active in the respective 
experimental situations but does itself not involve exci- 
tations of (p. The convection rolls might themselves be 
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smaller than the optical resolution and therefore remain 
unobserved. Then the onset of </? modulation would ap- 
pear to be the threshold of a primary instability of the 
homogeneous basic state, although the actual primary 
( "isotropic" ) threshold is at slightly lower voltages. 



F. Interaction with the pattern 

The equation of motion for the pattern amplitude itself 
is given by 



r(i9t+f7-V+ iv.xqc)A 



A. (10) 



All coefficients are real. The differential operators V^:, 
Vj,, Vz are used to indicate that only the immedi- 
ately following expression should be differentiated, not 
A. With the exceptions of VxdxA and Vx^xA, which are 
included for symmetry, only terms up to lowest nontriv- 
ial order have been included, assuming time, all lengths 
and all fields to scale independently. The l.h.s. is given 
by Galilean invariance, the form of the expression follow- 
ing from rotation invariance ]5^ . The term involving 
V z $r represents simply an additional contribution to the 
electric driving field. Symmetry would also allow a term 
Vy VyA, but tfic hydrodyuamlc equations do not gen- 
erate it. 

When Py — —Qilc it is possible to absorb the cor- 
responding term into the parenthesis after and to 
rewrite the complete expression in the "potential" form 
^y{dy — iqcip){dy — iqc(p)A. In this case the phase of the 
pattern does not drift in a weakly deformed ip field. In 
fact, (3y is close to this value (see Q, Appendix 0), 
which is largely due to the conservation of charge and 
momentum and the fiux-divergence form of the resulting 
expressions. 

The Landau coefficient g is of order unity in our nor- 
malization ll^. As in the conduction regime ji^, it 
is dominated by "geometric" effects, i.e., inhibition of 
the Carr-Helfrich mechanism for large Uz'. besides flexo- 
electric effects, about 99% of g come from contributions 
quadratic or cubic in Uz . This indicates that a breakdown 
of the weakly nonlinear expansion should be expected for 
A^nz^Oil). 

The nonlinear excitation of the soft modes by the con- 
vection pattern is less intuitive than their feedback onto 
the rolls discussed above. The principle of truncation for 
the contributions of A in Eq. ( [lO| ) is again to keep only 
terms of lowest nontrivial order, however, allowing for 
phase gradients of A without gradients of the modulus. 
As the relative scaling of ip vs. leng th sca les is left unde- 
termined at this stage (see Section IV D), the nonlinear 



term —A*iqcfA, which is given by rotation symmetry, 
must always come along with A*dyA. 

Since Eq. (40a) must have flux-divergence form, the 



largest contributions from A are O(V^A^), which is 
too small to be relevant. The coefficients Ir and /, in 



Eqs. (40t.40c) have the dimensions of a current den- 
sity and measure the strength of an alternating current 
j|^|2 = 2\ A\'^ {Ir cos ujt — lismujt) generated by the con- 
vection pattern. Important contributions to the coeffi- 
cients Ir, ... Jiz come from charge advection. 



The strongest contributions to T in Eq. (40d) are sim- 
ple potential effects. In MBBA, the most important one 
describes a relaxation of the bend of the director modula- 
tions by twist and is given by 4(fc22 — A:33) (see [ |33| , ^ and 
Appendix which is generally negative. The second 
most important contribution comes from the dielectric 
torques from applied and induced field on the director. 
In MBBA this contribution to F is positive and in ma- 
terials with large negative (e.g. ^ — e_L) it could com- 
pensate the elastic one and reverse the sign of F. Notice 
also that F is paricularly sensitive to flexoelectric effects 
(s. Table ||). A negative value of F is req uired for the 
occurrence of abnormal rolls (see Sec. VIA) and chevron 
patterns. 



The terms associated with Sx and Sz in Eqs. (40e 40g) 
represent internal stresses of the convection pattern. As 
for the coupling of A to in Eq. (p^, a correspond- 
ing term for the y direction or a term of the form 



Syxdylm{A* dx A} do not enter Eq. (40f), although they 
are allowed by symmetry. Sxx, Syy, . . . , Szx can be in- 
terpreted as surface tensions of the planes of equal phase 
of the convection pattern. 

The high number of soft modes and the rich, non- 
potential coupling almost certainly lead to spatio- 
temporally chaotic states already at onset, provided the 
spatial extensions of the sample are large enough. 



V. EFFECTIVELY 2D PATTERN DYNAMICS 
NEAR THRESHOLD 



Currently more interesting than the fully 3D chaotic 
state are, from the experimental point of view, the 
quasi 2D pattern dynamics in a restricted geometry near 
threshold [e' ^ (A/d)^]. The dynamics is here derived 
for a slight generalization of the usual setup: A magnetic 
field Hy might be applied along the y direction (i.e. in- 
plane, normal to the rubbing direction). With Hy slightly 
below the twist Freedericksz field XaHp = k22{'^/d)'^, the 
amplitude of the twist mode, which is known to be im- 
portant for the pattern dynamics in any case, becomes 
a slow variable and must be included explicitly in the 
2D formalism. The conventional setup is described by 
H„ = 0. 



8 



A. Derivation of the 2D description 

The appropriate boundary conditions for the soft 
modes at the enclosing electrodes are 



$0, V9, t/ = at 



±d/2. 



(11) 



Thus, all components of the electric potential and the 
velocity field are damped by the boundaries. 

Again the "center manifold" method is used (s. Ap- 
pendix P), now to reduce the 3D equations to 2D. The 
dynamically active part S of the state vector U is now 
given by the sum of 



A(x, y, z) — A'{x, y) cos(7rz/d), 
ip{x, y, z) = ip\x, y) cos{Trz/d), 
P{x,y,z) = P'{x,y), 



(12) 

(13) 
(14) 



with the amplitudes of the active modes A' , ip' , and P'. 

To be specific, associate A with the complex conju- 
gate of Eq. (|l^) and $o, ^r, ^i, ip, Vx,Vy, v^, and P 
with successive equations in the system (pO|), and define 
the scalar product (•!•) as usual as the equally weighted 
sum over z integrals over products of the two compo- 
nents. The projector onto the slow dynamics is con- 
structed from the bi-orthonormalized linear functionals 
(2/d) cos(7rz/d) Adz, (2/d) f^^^^^ cos{7t z / d) (p dz, 

and {1/d) J^^f^^Pdz. 

We proceed with a calculation of the excitations in th e 
fast subspace i? by a term-by-term solution of Eq. (B12). 



The truncation is chosen such that the distinguished limit 

(15) 



P\Hl- Hldt- s' 



p',A',d^,dy^e'^'\ 
e' = 0{X/df, and d ^ oo 



is correctly described. 

At linear order in the amplitudes and to linear order in 
dx,dy, the slaved part of the state vector contains only 
the contributions 



,(1) 



,(1) 







z2 








2 


/d2 








2 



d.P', 
d^P'. 



(16) 
(17) 



At order 0(|A'|^) and without any x oi y modulations 
there are excitations of the electric field 

$(2) + ) ^jz±M^ r \Aizr {\An^ dz', 

(18) 

and a contribution to the pressure field, orthogonal to 
the active pressure mode, 

p(^^^-{S.-2Se)-{\A\'^{\A\^)^), (19) 



where (|Ap)^ = \A'\'^/2 is the average of |Ap over z and 

Eoe_L{ui€^Ii + aj_Ir) , . 

OR := — ■ — ^ • 



By gradients of ^^,?\ P^'^\ and by direct contribu- 

tions at order 0{dx\A'\'^, dy\A'\'^) the mean flow 



.,(2) - 



16 



(5. + 5.-^.)^cos^(^) 
=(5.-5.)^cos^(!^) 



dx\A 



l\2 



TTZ\ 2dy\A'\^ 



(21) 
(22) 



is excited (there is no distinction between "singular" and 
"non-singular" mean flow in this approach). Excitations 
of $0 are of the order 0{e') and do, as all other remaining 
corrections, not contribute at leading order. 

Projection of the dynamics with the full state vector 
U = S + R onto the slow space yields the equations of 
motion for the pattern amplitude 



TdtA' = 



V'dy-^V' ] (23a) 



+ i—(iyVyip' + e' 

OTT 



4 

Se 



\A'? 
(23b) 

id^qcT %Se + (15 + 2^2)^^ - 95,^ , .,,2 



+ 



487r2 

id'^QcT 3 -I- 7r2 



127r2 

and the twist mode, 



A' 



lidtp' = {h^dl + k[^dl - Xa{Hl - Hi)) ^' 



2 

4d2 

r3 



Im <^ A' 
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3iqc I 



A' 



IT 

2f_ r 

37r3 



03 ^ 2a2 
m "4 



dxdyP' 



2a2{S_ 



E 



m 

and the pressure Poisson equation [|l[ 

d2 



dxdy\A'\^, 



0—12 
87r2 
47r2 



a4 



P' 



d^ SE + {l + 7ryS)Sx-S., , 

m 

d Se — Sz „2| Ai\2 



(23c) 
(23d) 



(24a) 
(24b) 

(24c) 



(24d) 

(25a) 
(25b) 
(25c) 
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B. Remarks on the 2D amplitude equations 



By the explicit representation of the couphng coeffi- 



cients in Eqs. (|23|-p5), the increa sing i mportan c e of the 
mean flow contributions in lines (23c, 23d, 24c, 24d ) as 
the separation d of the damping boundaries increases be- 
comes obvious. On the other hand, for small enough d, 
i.e. e' <ti (X/d)^, mean flow is negligible. With some 
rescaling (indicated by a caret ") dynamics are then de- 
scribed by the system 



fd{A = 



A, 



+ 2r Im{A*{dy ~i(f)A} 



(26a) 
(26b) 



This is the "normal form" for the dynamics of a pattern 
coupled to a soft mode with symmetry under the reflec- 
tions [x —X, (p —(p^ A A*] and [y —y, (p 
—if] or [y —y, Cp — > — (p, A ^^]- In the present 

case ci = C2 = I (^)^ ~ 0-96 are fixed by geometric 
constraints These values are quite close to the case 
with an underlying rotation symmetry ci = C2 = 1 [@, 
which turns out to be somewhat singular in its dynamical 

properties [^,|8). 

With i?^ = in line (24a), tp' is damped out and 
becomes one of many higher order corrections. Equa- 
tions (E3, Eq) with (^9 = are then sufficient, and with 
e' ^ (A/d)^^the simple Ginzburg-Landau Equation (|^) 
with 



(27) 



(the second contribution is numerically small, see Table |) 
is retained. 

When modulations of \A^ along x are strong, like in 
the chevron pattern, it is instructive to redefine the pres- 
sure P' P' — con st. X \A'\'^ such that its excitation by 
in line (25t) is canceled. It turns out that the re- 
maining excitation of P ' by dy\A'\'^ is proportional to Sx- 
The substitute for line (p3q) incorporating the correction 
in P' does then have the form 



zd2g,r(7r2 -Q) Se + S,- 



167r4 



and accounts for a mean flow with a nontrivial flow pro- 
file with zero z average. 



C. Pressure vs. singular mean flow 

As for any incompressible, Newtonian fluid, Eq. (|2 
is of the form 



In the present simple case, the matrix M is constant and, 
in the usual coordinates, diagonal. The inhomogeneity V 
(with dimensions of velocity) depends on time, cc, and y. 
Equation (^8|) can be solved formally by the transforma- 
tion 



dpP' = M.-^ izyidpG + V 



(29) 



which requires a fundamental solution G{x,y,t) satisfy- 
ing 

O^zxdp- dpP' 
= (z X dp)M.-^{z X dp)G + (z X dp)M-^V. (30) 

With some rearrangements, Eq. ( pO| ) has the same sim- 
ple structure as Eq. (|2|). The quantity G can be inter- 
preted as a stream function generating a certain com- 
ponent of the large-scale variations of the velocity field 
{vx,Vy) ~ (£ X dp)G, the "singular mean flow" (in early 
works |59| expressed in terms of the vertical vorticity 
-V^G). By using Eq. (|o|) instead of Eq. (|2|) and ehm- 
inating the pressure P' via Eqs. (|2^ ) also in the remain- 
ing equations [in our case ( p3| , p^ ] , a description fully in 
terms of the singular mean flow is obtained. 

In principle the three forms ( p8|p9|j30| ) are equivalent, 
although the last is often preferred in the literature (for 
an overview see both Refs. p],|60[|), perhaps because in 
some situations with high symmetry G is not excited, 
while P' is. When disregarding the effects of the addi- 
tional soft mode ip' , application of the less obvious but 
direct method of Kaiser and Pesch leads to the same 
result for the mean flow equation (pO) as the route de- 
scribed here. 

The flux V in Eq. (|9|) is determined by Eq. (H) only 
up to a transformation V ^ V + {zx df)Go{x, y, t) which 
implies a redefinition G ^ G — Gq. Hence (z x dp)G is 
not generally proportional to the z average of the large 
scale flow. This is not necessary for the formalism to 
work. In simple cases, like the present, the stream func- 
tion G has this property with the "natural" choice of V. 
To guarantee it in general, the method of Newell, Passot 
and Souh @ can be used to derive Eq. (El (Eq. (2.59) 



= dp-MdFP' -df-V. 



(28) 



in 1 60 1) directly under this additional constrain. 

Here, following Ref. |^] , the formulation as a mass con- 
servation equation (28) is used, since it derives naturally 
from the general formalism and is thus easily extended to 
three spatial dimension, to additional homogeneous soft 
modes, or to relax the assumption of incompressibility 
(which is not essential for the relevance of mean flow as 
is sometimes suggested). 

In principle, terms proportional to d'^cp', dyip' could 
also appear in the P' equation (^5|). It is a particular- 
ity of the system considered here that they do not. In 
the conduction regime of EC in cells with homeotropic 
boundary conditions it can be shown that there is 
such a, presumably small, excitation of P' by ip' propor- 
tional to the dynamic flexoelectric coefficient C,^ . 
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D. Variation of the boundary conditions 



VI. STABILITY OF IDEAL PATTERNS 



Of course, the reduction from 3D to 2D can also be 
carried out for other boundary conditions than (| 11). 
Some variants are of practical interest. 

An effect similar to applying a magnetic field along y 
can be obtained by a homeotropic (normal) anchoring 
of the director at the boundaries: For negative, not too 
small Ea, and with external electric fields as required for 
dielectric EC, a homeotropic director alignment is unsta- 
ble and instead the director becomes planarly oriented 
everywhere, except for small boundary layers of thick- 
ness (kii/eay^^E^^ « A. Hence, on the length scales 
^ A relevant for the amplitude formalism, the boundary 
layer vanishes and instead free boundary conditions for 
(p can be assumed. For this setup only some geometry 
factors have to be changed in Eqs. ([l0|-p5|). 

Another variant is a twisted cell, e.g. with = at 
z — —d/2 and tp = 7r/2 at z = d/2, as it was recently 
investigated experimentally by Bohatsch and Stannarius 
|l62| . Based on the reduced 3D description derived here, 
the linear theory for this configuration is developed in 
Ref. [lei. 



E. Higher order contributions 



It is worth noticing that, when only the limit of small 
e' is considered while keeping d fixed, there are several 
other nonlinear and higher order gradient terms besides 



those in lines (23c, 23d) which are formally of the same 
order of magnitude. A longer list containing more than 
fifty terms has been calculated numerically by Kaiser and 
Pesch ^l[ for the conduction regime of EC. Their form 
correctly predicts the stability of ideal roll patterns close 
to threshold but, since it does not separate mean flow and 
director effects, is unable to describe important effect s 
like the transition to abnormal rolls (see Sec. VIA). 
These limitations are partly overcome in the less system- 
atic but numerically surprisingly accurate description of 
Plant and Pesch @. 

In fact, since the lowest order mean flow effects enter- 
ing Eqs. ( p3[ - p5| ) all depend on gradients of \A'\'^, they 
do not contribute to long-wavelength instabilities of the 
band center (i.e. A = const.). For the calculation of 
the thresholds of long- wavelength instabilities of the pat- 
tern it may be useful to formally set up amplitude equa- 
tions including higher order mean flow. But for a sys- 
tematic quantitative description of general non-ideal pat- 
terns containing structures of size ~ Ae~^/^ only the trun- 
cation used here is justified. When higher order mean 
flow becomes relevant the 2D amplitude formalism is al- 
ready breaking down because, for example, the coherence 



length t.v£' 



'^1"^ becomes of the order of the sample thick- 



Rather than deriving stability bounds of ideal patterns 
A = a{z) exp{iqx + ipy) using a reduced 2D descrip- 
tion, it seems more appropriate to do the calculations 
directly based on the 3D equations, in particular when 
Hy ~ 0. This is easily seen from the fact that there 
is a (numerically small) manifest deviation of a{z) from 
the cos(7rz/(i) proflle at the threshold of all instabilities 
calculated below, indicating that the expansion for small 
A' is breaking down. Nevertheless, in order to obtain 
analytic estimates, Galerkin approximations will be used 
which correspond effectively to Eqs. ( p3[ - p5| ) and their 
extension to higher order contributions. For simplicity, 
the stability analysis shall here be performed only at the 
band center q,p = and only take homogeneous and 
long-wavelength instabilities with modulations along y 
into account. The latter restriction is justifled by ex- 
perimental observations and by the fact that, with this 
geometry, the advection of the patter by mean flow is 
particularly strong. It is then sufficient to consider only 
the interaction of if, Vx, and of the phase 6 given by 
A — a{z)expi6{y, z,t) \a — a {z) — 0(e' jg'Y^"^ be real 



and given by Eqs. (IC.lOb 40c)], which leads to the fol- 
lowing linear problem: 

71 a* = [(A:n + 2 be^^/e,) 9^ + 9^] ^ 
Pm. dt Vx = m {dy + dl) Vx + as dy dt (p 

+ Syy C? dy {OyO " Qc^f) (31b) 

+ S,, {a" 0) , 



2e)2/ 



rdtO = - TQcVx + Pydyif + iyd:^ 



(31c) 



Since a = at z = ±c?/2, the singular last term in 
Eq. ( ^lc| ) implies boundary conditions dzO = 0. The 
other boundary conditions are tp^Vx — at z — ±d/2. 



A. Homogeneous destabihzation 

First, consider homogeneous {dy — 0) destabilizations 
of the pattern. In this case (p decouples from Vx and 9. 
The destabihzation of (p is known as the abnormal-roll 
instability and was investigated in the dielectric regime 
in Ref. 1331. It was found that 



0-^ 9 fe33 

3d2g2r ■■ 



(32) 



which can be derived from Eqs. (23 2^), is typically a 



good approximation of the threshold. In MBBA e' 



AR 



2.8 X (A/(i)2 {ujtq > 8). The value o btai ned by numeri- 
cally calculating eigenmodes of Eq. (31a) directly is 3% 
lower than the value of Eq. (32). Below it will be shown 
that for MBBA the abnormal-roll instability is precede by 
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a long- wavelength modulation instability. Nevertheless, 
some phenomena associated with abnormal rolls might be 
observable around e = ear, e.g. the tendency of defects 
in the convection pattern to cluster along lines parallel 
to the rolls. 

For the discussion of homogenous perturbations of 
and 9 notice first the neutral mode associated with a 
translation of the pattern — > 6* + const.. This mode is 
best dealt with by decomposing 9 as 9 — 9{z,t) + 9(t) 
such that the z average {9}^ vanishes. By multiplying 
Eq. ( ^lc| ) by and integrating over z the dynamics of 
Q is obtained as 



dtB = -q. 



{a^Vx). 



(33) 



This describes the advection of the pattern by mean flow 
at large pattern amplitudes. 

At the threshold of instability one has dt9, dtv^ — 
and dtQ = const. {dtQ ^ implies an acceleration in- 
stability of pattern and liquid crystal). To calculate the 
threshold, eliminate 9 from Eq. ( |31h| ) by Eqs. (31c, 33|), 
obtaining the equation 



= mdlvx + 



(34) 



from which the critical mode can be determined numer- 
ically. When using the low-amplitude approximation 
— {e' / g') cos^ [it z / d) the critical mode is found to be 
antisymmetric in z (i.e. 9f0 = 0) at 



-drift — '''3-3 ^2 g 



[With the large amplitude approximation 



(35) 

= s/g, 

the first symmetric and antisymmetric mode both be- 
come unstable at the same e = {4:TT^rj2gS,'^) / {(P SzzTQc) ■] 
The antisymmetric excitation of the phase 9 involved 
in this instability is obviously inaccessible to a reduced 
2D description. Since Szz < in MBBA, the value of 
^drift ~ —100 X {ujTo)^^{X/d)'^ is negative and the insta- 
bility does not occur. However, since the elect ric contri- 
bution to Szz [the second term in formula ( A25)] is always 



positive and comparable in size with the hydrodynamic 
one, a positive Szz is thinkable for other materials. When 
Szz < 0, the mechanism leading to Eq. (|3^ ) is stabilizing 
- in particular for all perturbations of 9: The lamellae of 
EC are forced to align normal to the boundaries. This ex- 
plains why the experimental shadowgraph images, which 
average the pattern along z, remain quite sharp even for 
complicated pattern dynamics. 



B. Modulation instabilities 

The stability of ideal patterns with respect to pertur- 
bations in ip, Vx, and 9 modulated with small wave num- 
bers k along y was investigated numerically using the sys- 
tem (Sll). It was found that for MBBA the dominating 



destabilizing feedback loop is based on th e excitation of 
Vx by the term containing 5*^^ in Eq. ( |31tj ) and advection 
of the phase. There is a good Galerkin approximation for 
the numerical results. Using the low amplitude approx- 
imation for a^, Galerkin modes ip,Vx ~ cosjir z/d) and 
9^1, and projectors / cos{TTz/ d) ■ dz on Eqs. ( 31a , |311: ) 
and f co s^ (nz/d) ■ dz on Eq. (plq) to reduce the sys- 
tem (3^) to algebraic equations, the threshold for long- 
wavelength modulations instabilities is estimated to be 
at 



e — e' 



zz 



(Pqc (36) 
[512fc33^yyr - ri2T (256/3^ + 27n\,fy 



With Py = —qc^yy thc pareuthesls following F nearly 
vanish; the re main der is related to the small deviation 
of ci in Eq. (|26a| ) from unity. Numerically e'zz ~ 
9 X (a;ro)~^(A/d)2 (observe that e'^z and have differ- 
ent frequency dependence). In particular, at ujtq — 8, us- 
ing the second lowest Fourier approximation without flex- 
oelectric effects, Eq. yields e'^zT 0.792(A/d)2 while 
the numerical solution of system (31) gives a threshold 
at e' = 0.797(A/d)2. Using the Galerkin approximation 
it is easily verified that the instability is in fact of the 
long-wavelength (fc 0) type. 



VII. QUALITATIVE TRANSITIONS 



For the case Hp — 0, order of magnitude estimates 
shall be used to distinguish regions of qualitatively dif- 
ferent pattern dynamics in parameter space - above as 
well as below the stability bounds of ideal patterns. 

As mentioned before, simple Ginzburg Landau dynam- 
ics can be expected for small e until the last two lines in 
Eq. (HI) become relevant. Assuming dx ~ e'^^^^~^ and 
~ e' /g' , these terms have an effect of the magnitude 
of e'A' when, say, e' « [3 m g qc ^x A^] V[(^' " 6) {Se + 
Sx ~ Sz)Td^f « 2. X 10^ X (A/d)4(a;To)-2 (MBBA, 



> 10). Preliminary simulations of Eqs. (^J25) with 
= show that for higher e' defect cores (where Sa;!^!'!^ 
is large) are strongly deformed and lines along which the 
phase argyl' "jumps" are often generated and long living 
(rather long living phase jump lines are also observed ex- 
perimentally and in simulations of a similar model [|6ll ; 
it is not clear, though, whether these are due to lowest 
order mean flow effects). But the simulations also indi- 
cate that these lowest order mean flow effects, although 
they are formally dominating over the direct nonlinear 
saturation via the last term in line (23b), do not pre- 
vent the system form finally reaching a steady state with 
\^\^^e'lg'. 

Thus, assuming still \A!\^ « ^' 1 9' , there will 
be a further transition at e' « [3^/^77 (a4 -f 
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(MBBA, lutq > 8) where contributions from slaved exci- 
tations of the in-plane director (p by mean flow become 
relevant in the 2D pattern dynamics. Semiquantitatively 
these effects are described by Eq. ( p4| ) with Hy = 0, but 
a restriction of dynamics to a single linear mode of tp is 
then not justified. 

When, with increasing e, horizontal length scales be- 
come of the size of the sample thickness d [at e' « 
{T^ix/df « 0.2 X {X/df (MBBA), say] the 2D descrip- 
tion breaks down. Because then d is not the domi- 
nating length scale for the damping of mean flow any- 
more, the trend in the influence of mean flow on the 
smallest structures in the A field (e.g. defects) is re- 
versed. Now the structures themselves set the length 
scale. Assuming dx ~ e^^^^~^ and \A\^ « e/g, the con- 
tribution irVxA in Eq. ( p^ is large compared to eA up 
to e ^ {qcSxTix/mg? « 1- X 10^6 X (^^^)2 (MBBA, 
wTo ^ 10). For larger e the cores of defects are not 
affected by lowest order mean flow effects. Larger struc- 
tures, like the phase field of the pattern or variations in 
the defect density (e.g. in chevron patterns), may still be. 

Current experimental resolutions are of the order Ae — 
■ ■ ■ 10~^. The estimates above suggest that lowest 
order mean flow effects are best observed near the upper 
limit of the frequency range for the validity of the hy- 
drodynamic description used here, at wtq = O(IO^) (see 
Section [I C ) . With d « 5A they should be observable in 



the range of validity of the 2D description. 

According to the model for the chevron mecha- 
nism p5| , chevrons depend essentially on the abnormal- 
roll mechanism and can only form above the abnormal- 
roll instability bound [here given by Eq. (^)]. Thus, it is 
plausible to assume a e' ~ (A/d)^ threshold for chevron 
formation, in accordance with measurements presented 
in Ref. [ |2^ , where an approximate uj~^ frequency depen- 
dence is found. 

For very high e' many authors report the formation of 
disclination loops, which, being singularities in the direc- 
tor field, indicate already a breakdown of hydrodynamics. 
In thicker cells, the chevron pattern might decay along 
other routes before the 3D amplitude formalism breaks 
down at some e = 0{\). 

For convection in most quasi-2D systems, i.e. systems 
with d/X = 0(1), all these transitions, from the break- 
down of simple 2D Ginzburg-Landau dynamics to the 
breakdown of the amplitude formalism, do, in principle, 
collapse at e — 0(1). Only by specially designed experi- 
ments (as in [|66[) can these transitions be unfolded. 



expected for the unstructured basic state at a given volt- 
age, minus one. As for the Nusselt number in thermal 
convection, they are to first approximation proportional 

to \A\\ 

Several interesting questions can be addressed by mea- 
suring electric Nusselt numbers. First, it follows from the 
discussion above that, in dielectric EC, |Ap « e/g for 
thick enough cells even far inside the three-dimensionally 
chaotic range (the average of across the pattern is 
typically only weakly reduced in the presence of defects; 
see, e.g., Ref. [^). Measurements of the electric Nus- 
selt numbers therefore seem to be an effective method to 
test an essential feature of the theory. Second, Nusselt 
number measurements may also help to identify the qual- 



itative changes in the dynamics predicted in Section VI] 



in particular since they do, in contrast to optical meth- 
ods, not loose their sensitivity in thick cells or with small 
wavenumbers. Finally, the frequency dependence of the 
Landau coefficient g' , which enters the Nusselt numbers 
near threshold in a simple way, provides information on 
the strength of the dynamic flexoelectric effect (see Ap- 
pendix From the derivation of g' it is clear that 
boundary effects do not interfere in these measurements. 
To obtain the theoretical value for the Nusselt numbers, 
take the x and y average (symbol {■)xy) of Eq. (40b) plus 
i times Eq. (40c), which leads to 



{Ir+ll,)d.{\A\')^^ 

+ {Irz+ih.)d,{lin{A*dxA})^y. (37) 



Taking the boundary conditions (|^,pT|) into account, this 
implies, similar as for ^l?'' + in Eq. (O), 



^jz ■■= -(o-± + iujei_) ($r + i^i)xy 



+ {Irz^-lhz) {l^{A*dxA])^y^^ (38) 

ai z — ±(i/2. Obviously, jz is the complex amplitude 
of the average, pattern-induced electric current density 
through the sample. (The last term on the r.h.s. repre- 
sents a correction due to a global deviation of the average 
wavenumber from the critical one, and will be dropped 
below.) Thus, 



Mr 



(39a) 



VIII. ELECTRIC NUSSELT NUMBERS 

Recently "electric Nusselt numbers" Afr, Mi have been 
introduce by Gleeson, Gheorghiu, Plant [^ as the ra- 
tio of the in-phase or, respectively, out-of-phase compo- 
nents of the electric current to the corresponding values 



and 



(39b) 



At threshold, {d/de') (1^1^)^^^^ = {2g')~\ e.g. dMr/de' = 
0.31 in MBBA at ujtq = 8, dropping flexoelectric contri- 
butions. The value increases roughly ~ a; as frequency 
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increases. Typically, it seems to be larger than the cor- 
responding value for the conduction regime ^t\ . 



IX. CONCLUSION 

It has been shown how the 3D dielectric convection 
pattern interacts with various homogeneous soft modes, 
which are related to undamped hydrodynamic modes. 
The method to establish these relations is not unique, 
but the "center manifold" method seems to be favorable 
over the "order parameter" method. 

The reduction of the 3D pattern dynamics to a quasi- 
2D form in the layer geometry was derived analytically, 
thus establishing a description of the interaction of the 
pattern with the twist mode and the pressure field (or sin- 
gular mean flow). Scaling analysis suggests that the tran- 
sition from the simple, quasi-2D Ginzburg-Landau dy- 
namics to manifestly 3D dynamics in thick layers unfolds 
into several well distinguished steps, the first of which 
occurs already very close to threshold (e' = 0(A/(i)^). 



These characteristics should generally be expected for 3D 
patterns. 

Ideal, dielectric EC patterns are found to destabilize at 
some value e' ~ (a;To)~^(A/c?)^ for which an analytic ap- 
proximation in terms of material parameters is given. A 
particular nonlinear mechanism that stabilizes the phase 
of the pattern to be constant along z, thus giving the 
pattern a 2D appearance also at hi gher e', is identified in 
Sec. VI. As outlined in Sec. Vlll, measurements of the 



electric Nusselt numbers are suitable for quantitatively 
testing the theory, probing the dynamic fiexoelectric ef- 
fect in nematic liquid crystals independent of boundary 
effects, and investigating the route of the transition from 
simple 2D to fully 3D dynamics. 

It is my pleasure to thank L. Kramer, Y. Kuramoto, 
A. Lindner, W. Pesch, and E. Plant for valuable discus- 
sion, W. Decker and W. Pesch for providing the basic 
nematodynamic equations in a computer readable form, 
the Kyoto University for its hospitality and the Japan 
Foundation for the Promotion of Science (P98285) and 
the Ministry of Education, Science, Sports and Culture 
in Japan for their support. 
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Pm dt Vy 



= - [(^11 - 2baa) dl + {dl + )] $0 

- 2 be- {c7a/ea)dxdyip 

0= + Eodl [ey W - ((7,1 - Oa) ^r] 

- kg /eg + (1 - 3fe) ^0 ea/2 71] (e+ - 271 C"") $0 

+ [ (q-g/gg) (-fcll + fc22) 

+e_ (1 - 36) ( e+/27i - C^)] dy 8, ip 

+ (p-g/ eg) (^3 9^ Vx + ^2 ^z) 

+ iJo (/r dz 1^1' + /rx 9x Im{A* A} + 7^^ dz lm{A* 8^ A}) 
= -Eodl [ail $i + (e|| - Ca) a; 

- iJo (9^ + (a^ $i + e^ $,) 
- g; (1-26) (e+-2 7iC'')a^9,$o 

■ (-fell + fc22) - (2 6 e./e.) (e+ - 271 C"^)] 9. 8y 8, y 



7i8tip- 



Pm 8t = 



8tVz = 



+ ujdx (as dz Vx + a2 dx Vz) 

+ Eo (h dz \A\^ + hx 8x lui{A* 8z A} + 8z Im{A* 8x A]) 
: [/C33 8l + {hi + 26eVea) 8^ + ^2 5^ + XaH^] f 
+ (e+ - 271 - 26e_) 8x 8y $0 

- Q!3 8y Vx - Oil 8x Vy + {qcT/2)lm{A*{8y - iqcip)A} 
(ai +a4 + a5+ ae) 8^ Vx + m {9y + 8l) Vx 

+ (Q!2 + m) 8x {8y Vy + 8z Vz) 

+ 8x [-P - Ea e^ 8, + 2 C"" + 8^) $0] 
+ ag 8y 8t LP ^Sx8x\A^ + Sxx8xlvii{A* 8x A} 
+ Syy 8y\^{A* {8y " iQcp) A} + Szz8zlia{A* 8z A} 
(r]i 8l + ai8l + ai 8l/2) Vy 

+ 8y [{a2 + Vl) dx Vx + Q!4 8z V:,/2] 

+ 8y [-P-Eo€±8z^r + 2a2C^8l^o] +a28x8tp 
+ Sxy 8xlTa{A* {8y - iqc<p) A} 
[{a2+Vi) dl + ai (81/2 + 8^,)] Vz 

+ 8z [(03 + Oi2 + r]i) dx Vx + ai 8y Vy/2] 

+ 8z [~P - (fell - ^22 + 2 6e+ e_/ea) 8x 8y ip ] 

- [e+-2 6e+ -2(7i+Q2)C''] 5^ 9. $0 

- ^0 [(eii -ea)8l + e^8l + 2 e^ 8^] 

+ Sxz 8xlui{A*8z A} + Szx8zlm{A* 8x A} 
■■V -v 
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(Note: I^T^jiX is confused here, so we continue with a single 
column) 



APPENDIX A: RESULTS FOR LINEAR 
STABILITY AND COUPLING COEFFICIENTS 
IN 3D 

Some analytic and numeric results relating the 3D de- 
scription of the pattern dynamics to hydrodynamics are 
presented here, in particular analytic approximations for 
all cocfBcients entering the the results in Sections 0, [v| , 



and VII]. In the second lowest Fourier approximation 



the onset of dielectric EC is at 

4wCT|| {al - 71 771) 



2uja\\ (a2e|| +£0771) 



X \a2 .221 5 - fa ??l cr_L 



(Al) 



with a critical wavenumber 



2 2w (-a^ +7i?7i 
1c = — 



k33 m 

0,2,, 2, ^2 2eae|| r?i ay 

" " Q!2 A - Eg m Cr± 



2Xuj(T\\ (q2 e|| + 771) 

- (-a m O'-L 



(A2) 



where X :— £\\cra — eaf||- The accurate numerical result 
for vs. Lo is nearly a perfect straight line (see Fig. |^), 
which can probably be understood by means of the ap- 
proximation used in . Formula ( |AlD nicely estimates 
the offset of this line for intermediate wtq, but gives a 
different slope as wtq 00. The deviations at the lower 
end are an artifact occurring with all truncated Fourier 
approximations. For the situation is similar. 

In the second lowest Fourier approximation, the excita- 
tion of the hydrodynamic fields in the convection pattern 
near threshold (critical mode) is 

$ =|yl| sin(g,x-f argA) 

X ($„ -t- cos(2wi) - 2$, sin(2wi)) , ^ ' 
riz =2 |A| cos{qc X + a,ig A) {cos{uj t) + sm{uj t)) , (A4) 
/ =2|A| sin(g,x-f arg^) 

X {fc cos{ujt) - Js sin(wt)) , 

where / generates a velocity field v — {dxdz, dydz, —d^ — 
dy)f and the real constants fc, fs, ^u, ^c, are given 
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by 

(A6) 

4q!2/s9c = - ea^^c - + 4A;33gc 

eg^^e ~ ^11) + (^^ll'^ + ^11)] 

(A7) 

= (A8) 

9cO-|| 

With the help of the real constants $+, n+ 
which characterize the adjoint eigenvector 



^+ _ Eg {a2 e\\ + eg rji) Qc 



(AlO) 



$+ + i$+ = - X 

" ^ Aa2XEcu; 

{+E^g [-a2X + eam (2i 

-4 [(a^ -7i??i) w + fcssr/ig^] (2i e|| w + (7||)}, 

< = (A12) 

a2 

and the normalization factor 

4N=2{ql-nt)€aE^^ 

+ E, [{ea n+ - e\^ql) ($„ + + $,) (A13) 
-2e,^ ($+ + $+) _a„ (2 + $+-$+)], 

analytic results for most coupling coefficients entering the 
3D description can be obtained: 
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(A14) 



E. 



am (^2e[|ej_w^(T|| + 2e|W^(T^ + o-|0-±) - ^2 (2e|^^ + ^f) ^ } 
2e||A:33'7i ("2e|| - Ca??!) ^^9? + {(^2^ + '^aVi) ("2 - liVi) ^^o"!! + a2fc33'7i9cf^[| -'^ + ialE^.a\\X^ 

(A15) 

- 4 ($+ - $+ - 2 ($J + $+ $c + + flexo, 

27V^2^e_L£'c ($„ + $c + $,s) qc + ^a3Uj-2 [fc + fs) {ai + + + + {a^ / m) ll 
-Ae^LO ($+ - $+ - 2 a_L + $+ $c + *s) + flexo, 

(8/3) Ng = {2ql- Ant) E'l + 4 ql [uj (ag + 72) - 2 e„ n+ $e] 

+ (/e + Is) [(6 a2 + 4 aa) - A {ai + 72) 9^'] ~2E,q, (3 <i>^ + <i>^ - 2 $+) 
+ e,£;,g, [(5nJ - 2ql) ($„ + 2$,) - Auj (2$+ + $+)] + ^ojqle^ + 2$+$, - 2$+$,; 

+ 4 (7^ aa [2 + 2 + (2 flexo, 

4/^ + 4i/i (ieatj + cTa) [(4- 2i) Ec~{2-2i) q^ + - *s)] 

+ [(-3 + /e + (l-i) /s] + 2e||g^ [A ($„ + $c + i + » /. ($„ - $c - i , 

2 (ieati^ + CTa) 



(A16) 
(A17) 

(A18) 

(A19) 



AL 



Ail,, 



fa Ec 

472 (/c - /.) ql + 2ea£^c (*« + 3-, - - ($2 + 2$^ + 2$2) + Aq, - fcn) 
(2-2z) {ieaUj + aa){^u+i^c-<^s)-qleaE, [(-3 + z) /c + (l-i)/s] 
2 en g3 [/^ + + i + ifs ($. 



(A20) 



flexo, 



T^Ak22-Ak33-eaiE,/q,) ($„ + $,-$,) + £, ($2 + 2 + 2 $2) _ 2 (?c (/e-/,s)+ flexo. (A21) 



2 = [2 ql (ai + ^5 + 72) (/c - /s) - Ca E, ($„ + - + q, e,, ($2 + 2 $2 + 2 $2)] 
2S,=2a5ql {fc- fs) + eaEc [2 E, ~ q^ ($„ + $c -$,,)] , 
2^,, = (2ry2-a4) {fc - fs) ql + qc^i^ i<^l + 2^l + 2<i>l 



2S,,=2 (72 (/c-/s) g^ + ^cCi ($2+2$2 + 2<i>2 



2 'S't; 



(^2 + as) (/c - Is) ql - Ca (^-u + $c 



flexo, 
flexo, 

g,e|| ($2 +2 $2 + 2 $2)+ flexo. 



2^.. =2 (72 -ai) Uc-fs) ql + qce± {^l + 2 + 2 <S>1) + Aq, (^^33 - fcn) + flexo, 



(A22) 
(A23) 
(A24) 
(A25) 
(A26) 
(A27) 



Here "flexo" stands for flexoelectric corrections, which 
involve matrix inversions and are hard to express in a 
compact form. For similar reasons, no formulas are given 
for Kx, Sxx, Szxi Irz and Uz- The coefficients van- 
ishes for the second lowest Fourier approximation (and 
Pm = 0) and ax and (3y + qc^y have only flexoelectric con- 
tributions. The expression for has a different structure 
than the other formulas because it was not calculated 
with the "center manifold" formalism but by differentiat- 
ing the determinantal condition for stability of the basic 
state, which is more effective in this case. There are in- 
direct contributions from excitations of riz entering Sxz^ 



Ux and Irx [the bracket in Eq. (|A20|) , in Eq. (|A27|) they 
cancel favorably] . These have been calculated only in the 
lowest Fourier approximation (i.e. for the non-oscillatory 
part of 7I2), which introduces a small error. 

The quality of these results can be judged by compar- 
ing the values obtained for MBBA at wtq = 8 with the 
accurate numerical values in Table ^ (only the approx- 
imation for the indirect contribution from entering 
Sxzi hx and Irx is retained). In the combinations in 
which the results are expressed there, they are, for fixed 
ciJTo, independent of the electric conductivity and, with 
the exception of Kx, Kz, Irx , and some contributions in- 
volving ('^ , also largely independent of wtq for ujtq > 8. 
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Two exceptions, which should both be experimentally ac- 
cessible, shall be highlighted: the dynamic flexoelectric 
contribution to EcUx, which increases linear in lotq, and 
the contribution to g proportional to (^'^'^ which is the 
only one which increases ~ (wtq)^ as lutq —> oo. The 
result for g in the lowest Fourier approximation 



5 = 9.45 + 0.00252 - 0.00401 uj tq e+ C''^ 

-5.18X 10'"^ UJ tqC"^^ + 0.00102 oj^T^C'^^ (A28) 

illustrates the latter effect, although it is correct only in 
its o rder o f magnitude. The flexoelectric contributions in 
Eq. (A28) and also in Table || are expressed in terms of 
6+ := 61 + 63, 6_ := 61 —63, and (^'^ := 271^^, in units of 
lO^i^Cm-i (3.00 X lO^Myn^/^ in Gaussian units). Typ- 
ical values measured for ei and 63 are a few times that 
much (see the overviews in Refs. [|68| , p^ ). 

As a result of the approximate 7r/4 phase shift of the 
director oscillations, Kz is so small that the remaining 
finite viscous effect is comparable in size to the effect of 
finite mass density pm, which has been suppressed every- 
where else. 



APPENDIX B: "ORDER PARAMETER" 
"CENTER MANIFOLD" METHOD 



VS. 



Here, two general methods for obtaining amplitude 
equations are compared. It is shown that they give dif- 
ferent results in the presence of multiple homogeneous 
soft modes. In this work, the "center manifold" method 
is used to derive the reduced equations. In order to keep 
the formalism simple, it will be restricted to homoge- 
neous {q = 0) modes and their slow modulations. The 
inclusion of patterning soft modes {\q\ 7^ 0) is straight 
forward. 



1. Formal setting 

Let the state vector C/(r) describe the configuration 
of all relevant degrees of freedom (e.g. hydrodynamic 
fields) of the system in the (ideally) infinitely extended, 
D-dimensional r space. Assume the "microscopic equa- 
tions" to be of the form 

= F{U) = L{dt, V}U{r) + nonlinear terms, (Bl) 

where the linear operator L{dt,V) is polynomial in dt 
and V, acting on U{r) locally and translation invariant 
in space and time. There are several branches j of lin- 
ear modes Vj (q) which solve the generalized eigenvalue 
problem 



and, with some suitable scalar product (-I-) (which does 
not contain an integration over f), adjoint eigenstates 



W,iq)L{a,{q),iq)^0. 



(B3) 



For some branches j € K the growth rates Rc{(Tj(g)} 
vanish (or are small) at their maxima at g = 0, for the 
others (j ^ K) they are negatively large in the vicinity 
ofq = 0. 



2. The order parameter method 

Using the "order parameter" method physical states 
are characterized by weighted sums over slow eigenfunc- 
tions of L{dt, V), 

Uir) = J2 fu,{q)Vj{q)cxp{iq-r}dq + R{{uk}). (B4) 

The weights Uj{q) are interpreted as the Fourier trans- 
forms of the set of "amphtudes" used in the reduced de- 
scription. The range of integration 17 is a region around 
q = 0, large enough to include all significant contribu- 
tions from Uj{q) and small enough to exclude slow (pat- 
terning) modes at large wavenumbers. With this ansatz 
slaved contributions R{{uk}), which are fully in the fast 
eigenspace of L, come in only at nonlinear order. The 
linear dynamics for each amplitude is simply given by 



dtUjiq) = (J.j{q)Uj{q). 



(B5) 



An inverse Fourier transform yields the linear dynamics 
in physical space, which is usually simplified by truncat- 
ing the Taylor expansion of (Jj{q) in each component of q 
for small \q\, such that, in physical space, derivatives of 
Uj are obtained. 

However, the situation is different in the case of mul- 
tiple slow branches. Then aj{q) is typically non-analytic 
in the components of q (although it is analytic in |g|). As 
a generic example, consider the linear operator 



i(9t,V) 
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dxdy 



dt 



dxdy 

^dy- adt 



(B6) 



with a positive parameter a. There are two neutral modes 
at {q,p) := q = 0. It is easily seen that one of the two 
growth rates is of the form 



0-1 (9, P) 



P 



q^ 



q 



0(a), 



(B7) 



L{a,[q)MVM^O, 



(B2) 



i.e., non-analytic at q,p = 0. As a result, the correspond- 
ing amplitude equation in physical space 

7 

+ I Kif - r) Mi(f') d?r' + 0{a) + 0(^1,^3) (B8) 
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is nonlocal. In polar coordinates K{r) = 

— (3/7r)r~* cos4(/3. These conclusions do not require a 
to be small, because the additional terms 0{a) depend 
on a and they can cancel the non-analyticity calculated 
here at most at particular values of a. 

This transition from local basic equations to ampli- 
tude equations with algebraically decaying non-localities 
is counter-intuitive and misleading. This approach has 
the adv antage that in Fourier space the linear dynam- 
ics (B5) is simple. This is useful for calculations of pat- 
terns stability involving only a few Fourier modes. 

Finally, notice that a general method to reobtain lo- 
cal a mpl itude equations from the Fourier representa- 
tion (B5), e.g. by redefining the amplitudes, should not 
be expected. Equation ( p35| ) is general enough to include 
even non-local interactions in the basic equations, which 
certainly cannot lead to local amplitude equations. 



3. The center manifold method 

Alternatively, in the "center manifold" method only 
the slow modes at q = 0, Vj{0), are used for the charac- 
terization of the physical state 

U{f) = J2 f Uj{q)VjiO)eM^q-r)dq + Ri{uk}). (B9) 



The "slow subspace" spanned by the sum in Eq. (B9) can 
be extracted by the projection operator 



with 



dq 



dr' 



(BIO) 



{WM\ /(r'))F,(0)exp(zg. (r-P)), (Bll) 



where it is assumed without loss of generality that 
the states Wj{Q) and Vj{0) (j G K) entering P form 
a bi-orthonormal system. The "slaved" contributions 
R[{uk}) cover the remain ing subspace. 

The factor Vj (0) in Eq. ( B9) can be pulled out of the in- 
tegral, which is then simply the inverse Fourier transform 
of Uj{q) into physical space Uj{r). It is thus justified to 
define Uj{r) as (the local average of) the hydrodynamic 
variable {Wj\U{r)). In particular, if {Wj\U{'r)) is con- 
served, so is Uj(f). In such a local representation, the 
vicinity of the "center manifold" method to a multiple 
scale approximation would be more obvious. Similar sim- 
plifications would also be possible for the integrals below, 
but have been suppressed in order to ease the comparison 
with the "order parameter" method. 



The function R{{uk\) is defined by the perturbative 
solution of 



{l-P)F{U) = 0, 



(B12) 



where U is given by Eq. ( p39| ) and the Uk are small and 
vary slowly and smoothly in space and time but are oth- 
erwise arbitrary. 

At linear order in U , where Eq. ( |B12 ) reduces to 



(l-P)L(9t,V)C/ = 0, 



(B13) 



a form 



/ {q)Vj (0) exp(ig • r)dq (B14) 



with contributions rj{q) only in the vicinity of g = is 
sufficient. 

The amplitude equations are then given by 



PiL{du V)C/ = for each i G K 



(B15) 



where R is eli minated f rom U via Eq. (|B13| ). 

When Eqs. ( p313| ,B15) are satisfied with slowly varying 
Uk and small R{{uk\), this implies that U contains no 
fast eigenvectors of L(9(, V). Hence the resulting hnear 
dynamics for U is the same as the one obtained with the 
"order param eter" method, in particular R{{uk\) is then 
given by Eq. (|B14|) with 



J2 {Wk{q)\V,m-' {Wk{q)mO))m{q), 

(B16) 



k,l 



where I G K is running over all fast modes and j,k ^ K 
are running over all slow modes (for small \q\ the matrix 
{^k{q)\Vj{0)) in this expression is generally a perturbed 
unit matrix and readily inverted). 

To see that this method yields local dynamics for the 
amplitudes, split the linear operator, restricted to the 
subspace selected by Q := (1 — P), like 



QL{dt,V)Q = QL{0,0)Q + lidt,V) 



(B17) 



into a part Lq which is regular, and a term which is small 
for slow temporal and spatial variations of the operand 
and poly nomial in dt, V. Calling the sum on the r.s.h. of 
Eq. (|B9| ) S{{uk}), and supp ressi ng the arguments (such 
that U — S + R), equation ( p313 ) becomes 



QLR = QLQR= {Lo + l)R = -QLS, (B18) 
and is solved by expanding for small i.e. 

oo 

R=-Y.{-L^'l{dt,V))" L^^QL{dt,V)S. (B19) 



When ehminating R from Eq. ( B15| ) by Eq. (Bl£) and 
truncating at some power in the derivatives (i.e., for slow 
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enough variations), linear amplitude equations with local 
interactions are obtained. The exte nsion to the nonlinear 
level is straightforward (s. Section VA). 



It should be noticed that with this approach all modes 
in the kernel of L(0, 0) have to be treated as "soft modes" , 
some of which, e.g. those resulting from gauge symme- 
tries, may not actually have slowly relaxing modula- 
tions a ssocia ted with them (see e.g. th e pr essure mode 
in Sec. IV B ). For the simple example (B6) the "center 
manifold" method leads to amplitude equations identical 
to the basic equations. 

The reason for the difference between the two ap- 
proaches is that in the multidimensional kernel of i(0, 0) 
the choice of the basis vectors characterizing the slow 
modes is not unique. While they are fixed (with respect 
to the hydrodynamic variables in the "microscopic equa- 
tions") for the "center manifold" method, they point, de- 
pending on q, into arbitrary directions in the slow space 
for the "order parameter" method. For the same reason, 
there is no-near identity transformation mapping one rep- 
resentation onto the other. 
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MBBA 



lowest; 



4"^ lowest Fourier approximation 



Phase 5 



lowest 



qcPy 

Hz /to 
9 

Ir/EcUJe± 

Ii/EcUje± 
qcIrx/EcUje± 

qcIix/Ec<JJ£± 
qcIrz/EcLOe± 

qcIiz/EcUje±_ 

Se/Ec£± 

Se/qI 
Sx/ql 
Sz/ql 
S XX I qc 
^yy/ qc 

Szz I Qc 
Sxy / Qc 

Sxz/qc 
Szx/qc 

r 



7.52 X lO'' V^sm" 
1.53 X lO^'sm-^ 
1.01 
0.796 
1.39 + ... 

5.92 + ... 

-qUl + ■■■ 
-0.445 

9.37 + ... 

0.265 

0.0408 
-0.0229 + . . . 
-0.459 + ... 

0.7115 

0.2254 

0.0728 

16.6 
-2.23 

6.14 
-7.34 
8.32 + . . . 

-17.6 + ... 

9.86 + ... 
-13.7+ ... 

34.8 
-14.6 + ... 



7.74 X lO'-" V^sm 
1.79 X 10^ sm"- 
1.05 
1.06 

1.41 - 0.0089e^ + 0.0138e_e+ 

-0.0012e_C'^ 



0.0032e^- 



6.00 - 0.0031e^ + 0.0017e+C'^ 



+ 0.0005e+C -0.0004C 

/b2 



0.0014C' 



-g^Cy + 0.0022e_e+ + 0.0045e_C' 
0.0216 e+ +0.0439 C'^ 
-0.388 

6 X lO"'' - 3.13 X lO-^'m^ kg"^ p„ 



9.39 + 0.0056e3 



-0.0312 + 6.21 X lO-'^e^ 
-0.451 - 1.00 X 10' 



b2 



- 0.0452e+C'^+0.0755C' 
0.277 
0.0518 

3.12 X 10"^ e+ C'^ - 3.08 x 10"^ C' 
'e+C''^ +4.98 X lO-'^C'^ 



0.642 
0.250 
0.0851 

Below this line, all values are given in units of 10~^^N. 

17.1 
-1.96 

7.49 
-5.95 

7.07 + 0.0419 ei - 0.0561 e_ e+ + 0.0141 e\+ 

+0.0017 e- C'-^ - 0.0022 e+ C'-^ - 1. x lO^^'C' 
-15.7 + 0.0101 e- e+ - 0.0101 6+^+ 



+0.0050 e- C'-^ - 0.0052 e+ C'-^ - 1. 



8.66 - 0.0847 e_ 
-11.8 - 0.0201 e 



C'-^ + 0.0285 e+ C'-^ 



0.00500 c 



X 10~ 

lE'i 



b2 



-14.5 - 0.557ei +0.323e_e4 



0.0102 e+ C'- 
34.8 
0.0565ei 



1. X 10-" C""' 



-0.195e_C'^ +0.0665e+C'^ - O.OIOOC' 



b2 



4.01 X 10^ 



s m 



1.19 X lO-'sm" 
0.707 
1.89 
1.26 + ... 



5.59 



-0.690 

9.12 + ... 
0.293 
0.128 

0.590 + . . 

-0.502 + . 
0.992 
0.671 
0.162 

25.2 
-3.56 

18.2 
-16.4 

10.1 + .. 

-5.49 + . . 

10.2 + .. 
0.317 + .. 

88.8 
-32.1 + .. 



TABLE I. Coupling coefficients in the 3D dynamics of dielectric EC patterns at ljto ~ u!e±/a± = 8. The first column 
reproduces the analytic results Eqs. (A1-A21). The second and third column give accurate numerical results for MBBA and 
Phase 5 (Merck). For the units of the flexoelectric constants, see Appendix Ellipsis stand for suppressed flexoelectric 
contributions. 



"Using the parameter set "MBBA I" from Ref. Q. 
^As tabulated in Ref. @ for 30°C, e± = 5.22 eo. 



23 



2 




X {z+dl2) 

FIG. 1. Linear boundary layer calculated in the lowest Fourier approximation for MBBA. The temporally unmodulated 
component of the potential 3>„, the in phase and out of phase components of the director tilt nz,c, nz,s, and the in phase 
component of the velocity Vz,c are shown, normalization to unit slope at large distance (A <C {z + d/2) <C d) from the 
boundaries. 




24 



